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We analyze a quantum dot strongly coupled to the conducting leads via quantum point contacts - 
Fano regime of transport - and report a variety of resonant states which demonstrate the dominance 
of the interacting resonances in the scattering process in a low confining potential. There are resonant 
states similar to the eigenstates of the isolated dot, whose widths increase with increasing the 
coupling strength to the environment, and hybrid resonant states. The last ones are approximatively 
obtained as a linear combination of eigenstates with the same parity in the lateral direction, and 
the corresponding resonances show the phenomena of resonance trapping or level repulsion. The 
existence of the hybrid modes suggests that the open quantum dot behaves in the Fano regime like 
an artificial molecule. 
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I. INTRODUCTION 

In the past decades quantum dots were among the most studied systems in the solid state physic^^l. Proposed in 
the eighties as a system to minimize losses in the optical fibei'^, quantum dots have currently beco me a subject for 
fundamental research as artificial atoms^. Here the typical properties of an isolated natural atonpJ^ are qualitatively 
reproduced even in the presence of the interaction with the environmeniP. 

An artificial atom is a system created in a semiconductor heterostructure consisting of few electrons isolated from 
the environment by tunable barrierfP^. These non-infinite barriers allow for attaching conducting leads that open 
the quantum dot for transport whereas the properties of the isolated dot survive to a certain degree depending 
on the coupling strength. From the mathematical point of view, the quantum system admits on each transport 
direction a continuous energy spectrum with resonances instead of a set of discrete eigenenergies. In the case of very 
high barriers (weak coupling) the discrete energy levels of the isolated system turn into quasi-bound states. These 
are isolated resonances with an extremely small imaginary part, i. e., long life-time, whose real part can be well 
approximated by the eigenenergies of the isolated doiP. At weak coupling the physics of the transport phenomena is 
dominated by the electron-electron interaction that induces a shift of the resonant states'^ and the quantum dot follows 
the Coulomb blockade regim^3I!_ Decreasing the confinement barriers of t he dot the coupling with the conducting 
leads increases and both, the tunneling phenomena and the spin interactiorfiSini become more and more important 
relatively to the Coulomb interaction. In this intermediate regime the total transmission through the quantum dot 
shows broad and slightly asymmetric peaks, the so-called Kondo resonances,^ which are sensitive to the shape and 
the height of the confining barrier^^^. Upon further decreasing of the confining barriers the dot reaches the strong 
coupling regime. H ere the total transmission through the quantum dot shows asymmetric peaks and dips on a slowly 
varying backgrouncP^^. These peaks exhibit a Fano line shap^^ with a complex asymmetry parameter. They 
are narrower compared to the ones in the intermediate regim^. In the last years a number of studies were reported 
considering various specific aspects of transport in the strong coupling regime within non-interacting model d^ ' ^^'^ '' 
However ^ a s atisfactory theory providing a complete description of the scattering mechanisms in the low confinement 
potentiaP^at strong coupling does not exist. A particular difficulty, in this type of scattering potential is the high 
level density in the quantum system. As known from the particle physics, the methods used for describing quantum 
system s wit h a low level density, as the light atoms or nuclei, are not applicable for heavy nuclei with a high level 
densitjl22124| Por the mesoscopic physics, this means that the scattering problem for a quantum system in the strong 
coupling regime requires a different treatment compared to the scattering problem for a quantum dot in the Coulomb 
blockade regime. 

In the strong coupling regime, the electron scattering is profoundly affected by the quantum interference^25]_ The 
indistinguishabili ty of th e identical quantum particles leads to the interferenc^^I! between electrons and consequently 
to the Fano effecl!^^'^^'^. To explain this effect, often the existence of two interfering pathways or channels is invoked, 
one of which is resonant while the other is non-resonant. In the experiments in Ref. jl3j there are two spatially well 
defined interference paths consisting of the two arms of the Aharonov-Bohm ringSEni The arm in which a quantum 
dot is embedded defines the resonant path. In the experiments of Gores et aP there are no such clearly spatially 
separated interference paths, and the understanding of the Fano effect in this case is not straightforward. Clerk et 
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aP proposed as a nonresonant path a trajectory directly connecting the source and the drain contacts and as the 
resonant one a path passing through the dot via a resonant state and therefore spending a longer time in the dot. 
In this way, under certain conditions, the complex asymmetry parameter of a single resonance is associated with 
the dephasing time in the quantum system. In the frame of this model, the quite narrow and strong asymmetric 
("S-type" Fano) lines are found under the assumption that the quantum dot is coupled to two single- mode leads, 
but the slowly varying background is not explained. In Ref. |21j another model is discussed, in which the interfering 
paths are associated with the energy channels (subbands) of the leads: one of them contains a resonance or a group 
of overlapping resonances at the energy of the incoming electron while the second one contains only propagating 
states at this energy. As the result of the interference, asymmetric Fano lines with a complex parameter are obtained. 
In the presence of a scattering potential which couples the two channels it is shown that an interaction between 
resonances corresponding to different channels occurs, and this interaction exhibits dips in the total transmission for 
a favorable parity of the resonant states. As a second effect of the coupling between the scattering channels, the 
positions of the resonances corresponding to different channels are strongly modified in the complex energy plane. In 
the strong coupling regime the information about the scattering channels is actually not relevant for understanding 
the interaction between resonances. 

The above results confirm our earlier ideaP^to define the interfering paths using the resonances, i. e., the complex 
eigenvalues of the non-Hermitian Hamilton operator of the open systerrP^, instead of the quantum numbers of the 
lateral problem in the leads (scattering channel numbers). The resonances are also the singularities of the scattering 
matrix in the complex energy plane and, based on the decomposition of the S'-matrix in a resonant and a background 
ternPI, we have associated the interfering paths in a more formal way with these two terms. In the limit of a quasi 
one-dimensional model we have proven that the Fano function with a complex asymmetry parameter arises as the most 
general resonance line shape under the assumption that the background can be considered constant over the width of 
the resonance pole. The asymmetry parameter of the Fano line reflects the strength of the interaction between the 
considered resonance and the background which contains the contributions of all other resonances. These results were 
later confirmed in Refs. |19)22| . For decoupled scattering channels the Fano lines are only slightly asymmetric. The 
strong asymmetric ones, like those found by Gores et aF', imply the existence of many channels in the leads which 
can be coupled by dint of a nonseparable two-dimensional (2D) scattering potentia P^ - ^^ ' . As mentioned in Ref. ^2T] 
the interaction between channels changes the shape of the resonance dramatically. 

In this paper we develop a resonant scattering theory that takes properly into account the mentioned high level 
density in the quantum system as well as a strong coupling of the scattering channels in a nonseparable scattering 
potential. In view of Ref. [7] we assume that there exist direct trajectories which connect the source and the drain 
contacts, i. e., the potential energy in the region of the point contacts lies under the Fermi energy. According to the 
scanning electron microscope image of the systeirP, the quantum point contacts are very short and the leads are wide, 
allowing for a few subbands. The number of the conducting channels in the source and drain contacts is essential 
for the coupling mechanism of the quantum dot to the contacts. In the strong coupling regime and for a quantum 
system with a high level density, they limit the nu mber o f the eigenstates which couple to the continuum. The other 
eigenstates become consequently quasi-bound state^^HH!^ 

We believe that, for a deep understanding of the transmission through a quantum dot in the strong coupling regime, 
a resonant theory for two-dimensional systems is indispensable. In our opinion a resonant perturbation theory on 
the base of the Feshbach formalisrrPH is not sufficient for an accurate description of the strong coupling between 
the scattering channels. The resonances characterize the 2D scattering potential, and a direct solution of the 2D 
Schrodinger equ ation can not be avoided at least in the strong coupling regime. For this purpose we use here the 
R-matrix methocP^'^IHSI] and extend our scattering theory for ID systems without spherical symmetrjff^ to the case of 
2D systems. The R-matrix formalism is a very powerful method which allows for an efficient procedure to determine 
the resonances and for an exact decomposition of the scattering matrix into resonant and nonresonant contributions 
around each resonanc^^^. The second advantage of using the R-matrix formalism is that the scattering theory can be 
extended to describe the wave functions inside the scattering area. In this way the electron probability distribution 
density within the dot region can be analyzed, and the resonant states can be compared with the atomic orbitals. The 
coupling between the scattering channels leads to the occurrence of the hybrid resonant modes. Similar hybrid modes 
have also been evidenced in rectangular electromagnetic resonatorfP^l yielding a coupled mode with low radiation 
losses and a high Q-factor. As their atomic orbital counterparts, for example in _ff20- molecules, hybrid resonant 
modes arise in response to external perturbations of the isolated quantum system. While for atoms this perturbation 
consists of the molecular fields, for quantum dots the perturbation is the interaction with the conducting leads. 
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FIG. 1: (Color online) Potential energy in the quantum wire: constant potential energy in the source and drain contacts, Vi 
and V2, respectively, and position dependent potential energy in the scattering region {\x\ < dx, \y\ < dy). The quantum dot is 
isolated inside the quantum wire by the barrier (black area) with the height Vbo and the width db- The coupling between dot 
and contacts is set by the potential energy in the point contact regions, Vbi and Vb2- The constant potential energy felt by the 
electrons inside the dot is Vd- At the interface between different domains the potential energy varies linearly with position. 



II. THE MODEL 



We pr ovide a m odel for transport through a quantum dot that can be analyzed individually, like a single electron 
transistoj ^ l ^l^° l -^^ l The dot is embedded in a quite wide and infinitely long quantum wire, isolated inside a two- 
dimensional electron gas (2DEG) by infinite barriers, V{x, \y\ > dy) — )■ 00. Inside the wire the dot is defined by the 
barrier Vb(a;, y), which can have a general form (nonseparable and without any symmetry) like the black area in Fig. 
[1] At y = there are two quantum point contacts that ensure a strong coupling between the quantum dot and the 
rest of the wire, which plays the role of the source and drain contacts. The contacts are characterized by constant 
potentials Vs with s = 1 for the source and s = 2 for the drain. In the middle of the dot region the potential energy 
Vd is constant and can be varied continuously by a plunger gate. Further, we make the assumption that there are no 
bound-states in our system, i. e., min[Vi, V2] — imii[V{x, y)], \fx, \y\ < dy. 

The electronic wave functions are solutions of the two-dimensional Schrodinger equation 



2m* 



dy 



2 ; +y{x,y) 



ij{x,y) = E V(a;,y), 



(1) 



with the general nonseparable potential V{x^ y) in the dot region; E denotes here the kinetic energy of the electron in 
the plane of 2DEG and m* its effective mass. The Fermi energy^ for the 2D problem is fixed by the electron density 
Ns of 2DEG, Ep — nNsfi^ /gym* , where is the valley degeneracy factor. 

The electronic transport through a single quantum dot is essentially a scattering proces^^^'^ for which the potential 
energy has a spatial dependence only within a quite small region of the structure called scattering region (\x\ < d^^ 
\y\ < dy) and is constant outside it. As usual in the scattering theorjSSl this type of problem is solved using different 
methods for these two regions of the structure and the solutions are connected based on the continuity conditions of 
the wave function and its first derivative. 

Outside the scattering region the Schrodinger equation is exactly solvable and, as usual in the scattering theory, 
those solutions are superpositions of one incident and many scattered waves. 



4'i:\E;x,y) = 



0{N,{E)-n) 
2tt 



'Ssi exp [i kin [x + dx)] (t>n{y) 

00 

+ E Sln,in'iE)exp[~i kin' {x + d^)] (t>n'{.V), X < -d^ 

n' = l 

Ss2 exp [-i k2„ (x - d^)] 0„(?/) 

00 

+ E '5f„^2«'(£^) exp [i k2n' {x - d^)] 4>n'{y), X>dx 
n' = l 



(2) 



n > 1, s = 1,2, where denotes the step function, i. e., 6{t) = 1 for t < and 9{t) = for t < 0, and Sgs' the 
Kronecker delta symbol, i. e., Sss' = 1 for s = s' and Sss' = for s 7^ s'. The solutions ^ of the Schrodinger equation 
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are called scattering functions and the matrix S is the generalized scattering matri^^ or the wave transmission 
coefficients matrix. This is an infinite dimensional matrix, which connects incoming and outgoing components of the 
wave functions. T denotes here the transpose matrix. Due to the electron confinement in the infinite quantum well 
in the y-direction the functions 4>n{y) are given as 



4>n{y) = sin 



'dy 



Tin 

2dy^y^^y\ 



n > 1 (3) 



and the corresponding eigenenergies are 



2 



TT 



The quantum numbers n associated with the lateral problem define the energy channels for transport on each side of 
the scattering area, the so-called scattering channels. The wave vectors are defined for every channel (sn) as 



{E) = ko ^/iE- E^n-Vs)/uo, (5) 

where fco = TT/2dx and uq — h?kQ/2m* . In the case of the conducting or open channels, fc^n are positive real numbers, 
while for the non-conducting or closed channels they are given from the first branch of the complex square root 
function, ksn = Thus, the number of the conducting channels, Ns(E), s = 1,2, is a function of energy, and, 

for a fixed energy E, this is the largest value of n which satisfies the inequality E — E±n — > for a given value of 
s. The scattering functions exist only for the conducting channels. 

In the limit of a very low potential, i. e., V{x,y) — ^ 0, the scattering functions become plane wave corresponding 
to the free electrons. In the presence of a scattering potential with a non-separable character a plane wave incident 
onto the scattering area is reflected on every channel - open or closed for transport - of the same side of the system 
and transmitted on every channel - open or closed for transport - on the other side, the probability of each process 
being related to the elements of the generalized scattering matrix S. The 9 function in Eq. ([2| restricts the number 
of the elements with a physical meaning in S to Ni{E) + N2{E) columns for each energy E. 

For further determining the generalized scattering matrix S, the Schrodinger equation ([T]) should be also solved 
inside the scattering area. In this domain the potential landscape do es no t generally allow for analytical solutions and 
we have chosen to solve Eq. ^ by means of the R-matrix formalisnP32ll. Besides the extreme numerical efficiencjP^, 
this powerful method allows for a d irect comparison between the open quantum dot and its closed counterpart. In 
the frame of the R-matrix formalisnP^I^IHli]^ the scattering functions within the dot region, 

oo 

i,i^\E;x,y)^Y.^in(^)xi{x,y), (6) 
1=1 

\x\ < dx and \y\ < dy, are expressed in terms of the eigenfunctions xi corresponding to the quantum dot artificially 
closed by Neumann boundary conditions at the interfaces with the contacts. 



= 0, (7) 

x—:tdx 



dxi 
dx 

I > 1. Thus, the so-called Wigner-Eisenbud functions xii^^y) satisfy the same equation as ipn\x,y), Eq. ([l]), but 
with different boundary conditions in the transport direction: Since the scattering states 'ipn\x,y) satisfy scattering, 
i.e. energy dependent, boundary conditions derived from Eq. ([2| due to the continuity of the scattering functions 
at x = ±dx, the Wigner-Eisenbud function x;(a;, y) has to satisfy energy independent boundary conditions given 
by Eq. ([?]). The infinite potential outside the quantum wire requires Dirichlet boundary condition on the surfaces 
perpendicular to the transport direction for both functions, 'ipn\x,y — ±dy) = and xii^iV — ^dy) = 0. The 
potential energy for the Wigner-Eisenbud problem is given in Fig. |2][a) . As eigenfunctions of a Hermitian Hamilton 
operator the functions Xh ^ ^ 1: build a basis. The corresponding eigenenergies are denoted by Ei and are called 
Wigner-Eisenbud energies. They are real. 

The expansion coefficients a|^^ (E) are calculated using the Wigner Eisenbud eigenvalue problem and the boundary 
conditions satisfied by the scattering functions at the interface with the contacts. We have presented this method in 
detail in Ref. [S]- The coefficients a|^^(£;) are obtained as a function of S and the scattering functions within the 
dot region have the expression 

^{E; X, y) = -^@{E)[\ - {E)]K{E)R{E; x, y) (8) 
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FIG. 2; (Color online) (a) Closed quantum dot by means of Neumann boundary conditions at the dot-contact interfaces. 
This is the potential energy for the Wigner-Eisenbud problem, (b) Isolated counterpart of the considered open quantum dot 

(H = Hi = V62 ^ oo). 



with tpn\E;x,y) — ('^{E;x,y)j , n > 1, s = 1,2, and the R vector defined as 

V / sn 

RiE;x,y) 



Up Xiix,y)xi 



(9) 



1=1 



The vector xi in Eq. ^ is constructed using the Wigner-Eisenbud functions at a; = ztdj, and the eigenmodes 
corresponding to the lateral problem in the contacts, 



dy xiii-^ydx,y)(t>niy), 



n > I, s — 1,2. The wave vectors kgn define the diagonal matrix K, 

ko 



(10) 



(11) 



n,n' > 1, s,s' — 1,2. The matrix is also a diagonal one, defined as @sn,s'n'{E) = 9{Ns{E) — n) Sss' ^nn' , n > 1, 
s = l,2. 

Using further the continuity of the scattering functions on the surface of the scattering area one can derive a relation 
between the R-matrix 



R(£;) =uoJ2 



XI xF 
E — El 



and the generalized scattering matrix S, 



S{E) = [l - 2 (1 + iRK)"^l &{E) 



(12) 



(13) 



The equation ( 13 ) is the key relation for solving 2D scattering problems using only the eigenfunctions and the 
eigenenergies of the closed quantum dot [see Fig. [2ja)]. They contain the information about the scattering potential 
in the dot region and carry it over to the R-matrix. The matrix K characterizes the contacts and can be constructed 
using only the constant values of the potential in these regions. On the base of Eq. ( 13 1 the generalized scattering 
matrix S is calculated and further the scattering functions in each point of the system are determined using Eqs. ^ 
and The scattering theory together with the R-matrix formalism allows for a complete description of the open 
quantum dot and each physical parameter of the system can be further derived from the 5-niatrix. 

According to Eq. (12), R(-B) is an infinite-dimensional symmetrical real matrix and its expression allows for a 



very efficient numerical implementation for computing it. The big advantage of the R-matrix formalism is that, for a 
given potential landscape, only one eigenvalue problem with energy independent boundary conditions, i. e. Wigner- 
Eisenbud problem, has to be numerically solved and after that the generalized scattering matrix S can be constructed 
for each energy using Eq. (13). The computational costs are in this case minimal, but most important is that Eq. 



( 13 ) gives the explicit dependence of S on energy. This allows for an analysis of the scattering matrix and, after that. 



of the physical properties of the system, in terms of resonance energies. 
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The generalized scattering matrix S describes the scattering processes not only in the asymptotic region, but also 



inside the scattering area. But this matrix is neither symmetric nor unitary as can be seen from Eq. (13). Further, 
we define the current scattering matrix as 

The diagonal 0-matrix in the above expression ensure nonzero values only for the matrix elements of S that correspond 



to conducting channels for which the transmitted flux is nonzero. For simplicity, we have dropped in ( 14 1 the energy 
dependence of the matrixes and we will do this often henceforth. Using the i?-matrix representation of S, Eq. (13 1, 
the current scattering matrix becomes 

S = &[l-2{l + in)-^]&, (15) 

with the symmetrical infinite matrix 

"(^) = "oE|^ (16) 

1=1 ' 

and the column vector 

ai{E) = K'/'xi, (17) 



I > 1. According to Eq. (151 the current scattering matrix S is also symmetric, S = S . The restriction of iS-matrix 

to the conducting channels is the well known current transmission matrijP^EnUI]^ commonly used in the Landauer- 

Biittiker formalism. For a given energy E this is a {Ni + N2) x (iVi + N2) matrix which has to satisfy the unitarity 

condition SS^ = S = 1 according to the flux conservatioiP^. 

The elements of the current transmission matrix give directly the reflection and transmission probabilities through 

the quantum dot. For an electron incident from the contact s = 1, 2 on the channel n the probability to be transmitted 

2 , 2 

into the contact s' 7^ s on the channel n' is T„„/(£') — S2n'.in{E) ~ Sin',2n{E) ■ In the case of non-conducting 

(evanescent) channels these probabilities are zero. With these considerations, the total transmission through the dot, 
defined as 

Ni{E) N2iE) 

™= E E ^-'(^)' (18) 

n— 1 n' — l 

becomes 

TiE)=Tr[aiEW{E)], (19) 
where cr is the part of S which contains the transmission amplitudes, (Tnn'{E) — 82,1' .in{E), and cr^ its adjoint. 

III. RESONANCES 

The experimental analysis of a quantum system by coupling it to an electrical circuit has as a consequence the 
modification of its state. The physical interpretation of the measured quantities cannot be based solely on the 
properties of the isolated quantum system, but rather on the properties of the open system, i. e., the quantum 
system coupled to the contacts. Due to this coupling, the eigenstates become r eson ant states, some of them are long- 
lived resonances (called simply resonances) corresponding to quasi-bound state^^'^-' and the other ones are practically 
delocalized^ They can be found as states with a short life-time and it has to be elucidated if they influence 
significantly the p hysic al properties or not. In addition, the resonant states are eigenstates of the non- Hermi tian 
Hamilton operatoil^S^il of the open quantum system, and they are not orthogonal to each other anymor j^'^ l ^^ l In 
principle they can interact and their coupling may also influence the physical properties. 

From the mathematical point of view the resonances are assoc iate d with singularities of the current scattering 



matrix S. The representation of the iS-matrix in terms of $7, Eq. (151, allows for a very fast and efficient numerical 
procedure to determine its poles. When the quantum dot becomes open, the real eigenenergies of the closed system, 
migrate in the lower part of the complex energy plane, becoming resonant energie^^, Eqi = Eqi — iVi/2, I < 1. 
Based on this correspondence, we fix an energy E\ of the closed quantum dot and determine the resonance energy 
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Eqx as a solution of the equation det[l + in{E)] = 0. The matrix il, Eq. (16), is spht into a A-dependent part and a 



rest rix which should be a slowly varying energy function at least around an isolated resonance, 

niE)=uo4^+f^xiE)- (20) 

As shown in Appendix A, this decomposition of H around the resonance A leads to an expression of the S matrix in 
which the resonant and the background parts are separated 

SiE) = 2»"o //^^^_®-^ +^A(g), (21) 

where 

ME)^{l + inx)-'dx (22) 
is an infinite column vector that characterizes the resonance A, 

£x{E) ^ ~idx ■ (23) 
is a complex function which assures the analyticity of the current scattering matrix for every real energy and 

Sx{E)^&[l-2{l + inx)-^]& (24) 
is the background matrix. We have already proposed in Ref. |16j a decomposition of the scattering matrix similar to 



Eq. (21), but for an effective one-dimensional scattering system without channel mixing. There, the fJ-matrix is a 
2x2 one and the inversion of 1 + ifl reduces to a simple algebraic calculation. In the presence of channel coupling 
the inversion of an infinite matrix was a real mathematical challenge (see Appendix [A|). 

Based on the expression ( pTj ) of the scattering matrix, a resonant theory of transport through open quantum systems 



can be developed. Eq. (21) allows for the calculation of the resonance energies and for the analysis of each resonant 
contribution to the conductance. A similar decomposition of the transmission coefhcient in a resonant term and a 
background is also proposed in Ref. |22j . but in that case the two contributions can be evaluated provided that the 
eigenvalue proble m fo r the effective Hamiltonian of the open system is already solved. 



Based on Eq. (21), the position of the resonance Eqx — Eqx — iT/2 in the complex energy plane is given as a 



solution of the equation 

E-Ex-£x{E) = Q, (25) 
which can b e so lved numerically very fast using an iterative procedure starting with E — Ex- The complex function 



£x{E), Eq. (23 1, contains contributions from all Wigner-Eisenbud energies and from all scattering channels, i. e., all 
matrix elements of K. Thus, the resonance energy Eqx can totally differ from Ex and only in the case of a very low 
coupling of the dot to the contacts. Ex can properly approximate the real part of the resonance energy Eqx- With each 
resonance one can associate a resonance domain, which is a circle of radius Tx around Eqx in the complex energy plane. 
The resonance energies for the quantum dot shown in Fig. [T]are plotted in Fig. [3] together with the Wigner-Eisenbud 
energies. The geometrical parameters used for the numerical calculations were taken from the electron micrograph 
of the single electron transistor (SET) analyzed in Ref. [6]: 2dx = 2dy = 175 nm, db — 35 nm so that the electrons 
are confined within a domain of about 100 nm in diameter and the point contact regions are about 35 nm x 35 nm. 
The density of 2DEG is Ns = 8.1 x 10^^ cm~^ and the Fermi energy Ep — 29,6 meV. The confining barrier was 
considered VbQ = 100 meV and in the point contact regions it was taken Vbi = Vb2 — 2.5 meV; In the source and drain 
contacts Vi ~ V2 = 0. At each interface between two domains the potential energy varies linearly within a distance of 
10 nm. For the numerical calculation we fixed the number of the scattering channels to Ni = N2 = Np = 12, where 
Np is the number of the conducting channels at the Fermi energy. The analyzed scattering potential is not attractive 
and, in turn, it is not expected that the evanescent channels play an important rold2ni_ Numerically, the inclusion of 
the evanescent channels has not produced significant variations of the conductance. In Ref. [31j a detailed discussion 
is presented about the scattering potentials that allow for evanescent modes and about the influence of these modes 
on the total transmission through an open quantum system. 

As can be seen from Fig. [3] the open quantum dot strongly coupled to the source and drain contacts supports 
resonances with different widths, from very narrow, generally associated with m odes localized within the dot region, 
to very wide. This phenomenon is known in the literature as resonance trappin^^^- only certain states of an open 
quantum system with overlapping resonances couple with the environment and their widths increase with increasing 
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FIG. 3: (Color online) Resonance energies Eo\ = Eo\ — iV\/2 (red filled symbols) of the open quantum dot given in Fig. [T] 
with Vd ~ 0.0255124 eV and real eigenenergies Ex (black empty symbols) of the corresponding closed dot. 



the strength of the coupling, while the other ones are more or less decoupled from the continuun j^^ l ^^ l As illustrated in 
Fig. [3j for the quantum dot considered here there exist a few resonances with a long lifetime, but the majority of the 
resonant states couple to the contacts. This process is controlled by the number of the conducting channels according 
to Refs. |23|24j : In the strong couphng regime Ni{E) + N2{E) resonant states couple to the environment becoming 
quasi-delocalized. As we have shown in Sec. |ITj this number is energy dependent and increases with increasing E. 
This result is physically correct because the poles of the scattering matrix having real part much higher than the 
scattering potential, E ^ max[I^(a;, y)], have also a large imaginary part irrespective of the potential landscape. 

The decomposition (21) of the scattering matrix S in a re sonant and a background term is especially relevant for 
energies inside a resonance domain. According to Eq. (21 1 all matrix elements of S and, in turn, all transmission 
coefficients T„„' between the scattering channels have a similar dependence on energy around a resonance. In Fig. [4] 
we plot the transmission between the channel (11) and the channels (2n), n < Np for energies around a fixed isolated 
resonance. In the case of a symmetrical system the parity plays an important role. For the odd quantum number n 
the function (j^niy) has the same symmetry as (j)i{y), and the transmission coefficient Ti„ has a maximum around the 
resonance. If the parity is not conserved the transmission is forbidden, i. e.. Tin — 0, n = 2,4, .... The plots in Fig. 
|4] confirm the similar energy dependence of the transmission coefficients and we can conclude that a resonance can 
be completely characterized by the sum of these coefficients, i. e., by the total transmission. In Ref. [20l a similar 
idea was proposed and a global Fano asymmetry parameter was defined as a linear combination of the parameters 
corresponding to different scattering channels. 

The plot of the transmission coefficients. Fig. |4j shows a strong coupling between the scattering channels in the 
Fano regime of transport. The two quantum point contacts, specific for the SET geometrjl^, control the strength 
of the coupling with the rest of the quantum wire and confer the scattering potential its nonseparable character 
responsible for the channel mixing. In this case, a resonance cannot be associated anymore with a single scattering 
channel as proposed in the models based on the Feshbach formalism in Refs. |21|40j . The resonance perturbation 
theory'^ used there to describe the coupling between the scattering channels can have limitations for large coupling 
strength and becomes certainly very laborious for a system with many conducting channels. In our model the 2D 
Schrodinger equation is directly solved combining the scattering theory with the R-matrix formalism and this method 
can be used to describe each coupling regime. 



IV. CONDUCTANCE THROUGH OPEN QUANTUM DOTS 



The most common method to an alyze experimentally a quantum dot is to measure its conductance. In the limits 
of the Landauer-Biittiker formalisnF^^ and for very low temperatures, the linear conductance is given as the total 
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FIG. 4: Transmission as a function of energy around the resonance energy Eqx — 0.026537 eV for the quantum dot described 
in Fig. [l]and Vd = 0.0255124 eV. 



transmission through the dot at the Fermi energy, 

GiVa)^^TiEF-Va), (26) 
n 

for different values of the potential energy in the dot region. Each variation of Vd changes the scattering potential 
and, in turn, the total transmission. 

In Figs. [5f^a),[6][a),[7][a) the conductance is plotted as a function of Ep — Vd for the quantum dot presented in Fig. [l] 
with the parameter given in Sec. |III[ The conductance shows peaks with line shapes from symmetric Breit-Wigner up 
to strong asymmetric ones and even dips or antiresonances. These maxima and minima are usually associated with 
resonances. Some peaks in the conductance reach values greater than 1 and that means that at least two resonances 
interplay to determine the line shape. 

We analyze further in detail, in terms of resonances, each type of peak and dip of the conductance. For this 
purpose we need a functional dependence of the total transmission on Vd, at least an approximation, around the 
peak maximum Vq. In the case of a dip in the conductance, Vq denotes the position of the minimum. The R-matrix 
formalism used for solving the scattering problem allows, in a sense, for a very intuitive approach of T{Ep,Vd)- A 
small variation 5V = — Vq of the potential energy felt by the electron in the dot region can be approximately 
seen as a shift of the potential energy in the whole scattering area. In turn, the Wigner-Eisenbud energies are shifted 



with 5V and the Wigner-Eisenbud functions remain unchanged. For the R-matrix, Eq. (12 1, we can then write 
R(i?; Vb + 5V) ~ R(i? — 5V\ Vq)- This approximation is also valid for the total transmission 

T{Ef;Vq + 5V)^T{Ef-5V-Vo), (27) 

because the wave vector ksn is a slowly varying energy function. A detailed discussion about this approach is given 
in Appendix A, Ref. |16) for the open quantum dot without channel mixing. 
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FIG. 5: (Color online) (a): Conductance as a function of the potential energy in the dot region. The maps represent the 
electron probability distribution density inside the dot, \^n\EF;x,y)\^ for Vd = Vj}"""'""' (vertical dashed lines) for which 
j^(n^, riy) ^ rpj^g incident scattering channel is n = 1 for the odd modes in the lateral direction and n = 2 for the even ones, 
(b): The eigenenergies E„^^„y and the maps of the eigenstates, \ipn^,ny{x,y)\'^ , for the isolated dot. Bright corresponds to high 
values and dark corresponds to low values. 
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FIG. 6: (Color online) (a): Conductance as a function of the potential energy in the dot region. The maps represent the 
electron probability distribution density inside the dot, \^n\EF;x,y)\^ for Vd = yj"^'"!') (vertical dashed lines) for which 
j^{n^,ny) ^ rpj^g incident scattering channel is n = 1 for the odd modes in the lateral direction and n — 2 for the even ones, 
(b): The eigenenergies E„^^„y and the maps of the eigenstates, \4>n^,ny{x,y)\'^ , for the isolated dot. Bright corresponds to high 
values and dark corresponds to low values. 
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FIG. 7: (Color online) (a): Conductance as a function of the potential energy in the dot region. The maps represent the 



electron probability distribution density inside the dot, \'p^\EF\x,y)\ for Vd = Vq ^' " (vertical dashed lines) for which 
JSq"""'""' ~ Ef\ The incident scattering channel is n = 1 for the odd modes in the lateral direction and n — 2 for the even ones. 



(b): The eigenenergies En^,ny and the maps of the eigenstates, j-i/'n 
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Based on the relation (27), each peak in conductance can be associated with one or more resonances. We consider 
first an isolated resonance with the complex energy Eqx = Eq\ — iT\/2. The total transmission shows a peak around 
Eq\ that can be also seen in conductance if Eqx matches the Fermi energy, Eq\ = Ep. Therefore, a maximum in 
conductance at Vd = Vqa corresponds to a resonance A and the quantity Ep — Vqa gives the position of the resonance 
energy with respect to Vd- In this way the resonance energies can be directly compared with the eigenenergies of 
the isolated quantum dot. In view of the experiments presented in Ref. 0, this is a square dot with the dimension 
2d X 2d, d — d^ — dh, confined by a hard wall potential as depicted in Fig. [2Fb). Its eigenenergies 



TT 



2 



En^^n, ^^d+^ [Yd) + '^y)' ''^ - ^' 

are plotted in Figs, [sjb), [6j^b), [7]^b). The positions of the resonance energies are indicated by dashed lines in Figs. 
[5]ja),[6ja),[7]ja). As can be seen from these plots, the open character of the quantum system determines a shift of the 
eigenenergies in the complex energy plane, not only on the imaginary axis but also on the real axis. Due to the two 
quantum point contacts, which couple the quantum dot to the source and drain, the symmetry of the square dot is 
broken and the level degeneracy for rix, 7^ riy is lifted. 

A deep understanding of the transport properties through the open quantum dot requires a detailed analysis of 
the electron probability distribution density within the dot region, and the comparison of the resonance energies 
of the open dot with the eigenenergies of the isolated dolP3. In the upper part of Figs. |5|a), |6];a), and[7];a) the 
functions a;, y)^, n = 1 or n = 2, are given for x and y inside the scattering area and for Vd corresponding 

to the maxima and minima in the conductance. These functions are called resonant states or resonant modes. For 
comparison the eigenstates |Vi„^_„^ (x, of the isolated dot [see Fig. [2|b)] are presented in Figs. [5|b),|6|b), and 
[7]jb), where 



1 



d 



2d 



iJn^.nAx^y) = -sm -^(x + d) sin -^(y + d) , (29) 



2d 



n-x > i and Uy > 1. The function \ipn.^.ny{x, y)|^ has maxima in the x-direction and riy maxima in the y-direction. 
All modes {n,j;,ny) that we know from the isolated dot are also found for the open dot. Some of them are strongly 
modified due to the coupling with the contacts, but there are also modes that do not change much. Based on the 
similarities of the scattering functions at the resonance energy, \ipn\Ep; x,y)\'^ for Vd — Vqa for which Eqx — Ep, to 
the eigenfunctions \ipn^.ny{x,y)\'^ of the isolated dot, we associate further a pair of quantum numbers {nx,ny) with 
each resonance A, and the resonance energies Eox will be further on denoted by = _ ir("-"«)/2. 

The potential energy in the dot region Vqa, for which E^"'""^^ matches the Fermi energy, will be denoted by yj'"^'"«\ 
In this way the resonances are classified using a very intuitive criterion. 



A. Peaks associated with isolated resonances 



First we analyze the slight asymmetric conductance peaks associated with an isolated resonance denoted by A or 
by {nx,ny). In the energy domain of this resonance the scattering matrix S is given as a sum of a resonant term and 
a background, Eq. (21). Based on this relation and on the definition ( jTo] ), the total transmission can be similarly 
decomposed. According to the relation ( [27| and for small variation 5V of the potential energy around Vqa (for which 
Eqx — Ep), the conductance, Eq. ( [26| , follows the energy dependence of the transmission and becomes 

(30) 



GiVox + SV) ~ GresiEp - SV; Vox) + Gi,g{Ep - 6V; Vox). 



The resonant contribution to the conductance is an energy dependent function defined for each value of the potential 
energy in the dot region as 



Gres{E; Vd) 



2e2 



Tox{E) 



2i 1 


2 


1 


2' 


E-Ex-£x{E) qx{E) 




qx{E) 





with 



Tox{E) = 

and the energy-dependent Fano asymmetry parameter'^ 

1 

qx{E) " Jh.' 





bo 




PlX 




/32A 



-/3{a'^a/32*a, 



(31) 



(32) 



(33) 
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where (/^sa)™ = {l3\)sm s = 1,2 and {(Tx)nn' = S2n,in', n^n' > 1; The symbol * denotes the complex conjugate. The 
background contribution to the conductance is given as 

G,,{E) = ^Tr[cT,iE)aliE)]. (34) 



The functions Gres and Gbg are obtained from the expression (21 ) of the scattering matrix without any approximation. 

The first contribution to the conductance, G^es, contains a resonant term singular at E = Eqx and a term l/q\ that 
describes the coupling of the resonance A, characterized by the vector f3\, to the other resonances, characterized by 
the background matrix (T\. The function Gres yields always a peak mainly localized in the resonance domain. Due to 
the coupling of the considered resonance with the other ones, this peak cannot have in principle a Breit-Wigner line 
shape, even in the case of a narrow and isolated resonance; The lowest approximation for a resonant peak is a Fano 
line shape with a complex asymmetry parameter obtained for qx{E) ~ constant. The two terms add coherently to 
the conductance and it is usual to call Gres the coherent parll^^. The second contribution to the conductance is the 
noncoherent partP^ given only by the background matrix (T\{E). In the case of an isolated resonance A, it is expected 
that Ghg is almost constant inside the resonance domain. 

The conductance curve given in Figs, [sja), [6]Ja) and [Wa) shows two peaks that can be associated with isolated 
resonances. They correspond to the resonances (1,1) and (2,2) as follows from the analysis of the electron probability 
distribution density in the dot region [first and third maps in Fig. [sja)]. For these two peaks the resonant and the 
background contributions to the conductance are plotted in Fig. 8. As expected, the resonant part is given by a 
slight asymmetric Fano line and the background is almost constant. But, unexpected is the fact that the two peaks 
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FIG. 8: Conductance peaks associated with isolated resonances: {a){nx,ny) — (1,1); {h){nx,ny) = (2,2). Upper part: Resonant 
part of the conductance Gres; Middle part: Background part of the conductance Gbg- Lower part: Poles and position of the 
Fermi level in the complex energy plane. The potential energy in the dot region is constant, Vd = Vj}""^'""' 



are quite wide compared to the other ones in the conductance curve. The only possible explanation is related to 
the presence of the quantum point contacts, which modify dramatically the scattering process and the picture which 
we have from the effective ID scattering problem is no longer valid. In the case of the quantum dot studied here, 
the coupling between the scattering channels dominates the transmission through the dot and the scattering problem 
cannot be anymore reduced to a series of ID problems. In turn, in the presence of the channel mixing the resonance 
widths do not increase monotonically with the energy. 

As can be seen in Fig. |8j near the main resonances (1, 1) and (2,2), there exist other ones denoted by "a". They 
are broader, i. e., larger imaginary part, and are associated with modes localized mainly in the region of the two 
quantum point contacts, as shown in Fig. [9j The resonances of type "a" do not influence directly the transmission 
through the open quantum dot, but they play a decisive role in the coupling process of an eigenstate to the continuum 
of states in the contacts. There are many modes localized in the point contact regions with the resonant energies 
around the Fermi energy, but only for a favorable symmetry they can intermediate a coupling between the quantum 
system and the source and drain contacts. The probability distribution densities given in Fig. [9] shows evidently a 
coupling of the modes of type "a" with the resonant modes (1, 1) and (2,2). In this case, one can speak about an 
interaction between the two types of resonances. In the next section we will study this phenomenon for resonances 
localized within the dot region, which are very close in energy and have the same symmetry in the lateral direction. 
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FIG. 9: (Color online) Electron probability distribution densities: (a) (s, and (b) \il)'"^\x,y)\'^ for the mode associated 
with the resonance "a" in Fig. [s^a); (c) |i/)2^^(a;, and (d) \'tj}^\x,y)\'^ for the mode associated with the resonance "a" in 
Fig. ^h). 



B. Peaks cissociated with overlapping resonances 

Even in the case of a simple dot geometry, there exist only a few isolated resonances. The other peaks in the 
conductance with strong asymmetric line shapes or maximum values about 2 are typical for the scattering processes 
dominated by two or more resonances, whose resonance domains cross, i. e., overlapping resonances. How strong 
the overlapping resonances interact is determined by their relative position in energjED and by their symmetry in the 
lateral directioiPi] (perpendicular to the transport direction) . 

Let assume that there exists a second resonance A' = (n^,n^) around the Fermi energy, i. e., in the vicinity of the 
first resonance A = (nx^ny), and this is a broader one, > F^. The presence of the second resonance leads to a 
strong variation of the background term S\ with the energy around Ep. The expression of S\, Eq. (24), similar to S, 



Eq. (21 1, allows for a further decomposition of this term in a second resonant term and a new background. Following 
the method described in Sec. IIIIl we can write 



SxiE) = 2mo 



-S',{E), 



(35) 



where £'x{E), and can be obtained from (3\, SxiE), and Sx, Eqs. (22), (23 1, and (24), respectively, by 

;^ 

replacing dx by dx' and fix by $1^ 



— uq e'-eK • Thus, the background contribution of the first resonance to the 
•ibutions, 

(36) 



conductance, Eq. (34), becomes a sum of two contributions, 

GbgiEp', Vd) = G'^esi-Ep] Vd) + G[g{Ep; Vd 

a resonant one, 



G'resiE; Vd) 



2e^ 



Ux{E) 



2i 1 
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E-Ey-£'x{E) q'xiE) 









with a similar energy dependence as GresiE', Vd) and a second background, 

G',g[E:yd) = %-TT{a',^{E)a'l{E)\ 



(37) 



(38) 



with (cr^)„„; = In' I ^t"^ ^ Ij slowly varying with the energy if a third resonance does not exist around Ep. The 
energy-dependent Fano asymmetry parameter q'xi^) ^'i- (37), associated with the resonance A', has the expression 



(39) 



OA 



where Tq^{E) = l3[^ P2X and {(i'gx)n = {P'x)sn, s = 1,2, n > 1. The expression (|36| is also exact and we have 

only rearranged the terms in order to put directly in evidence the contributions of each resonance to the conductance. 
The second background term, G'^^g, gives the possibility of a further decomposition in a third resonant term and a 
new background in the case of three interacting resonances around the Fermi energy. 
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For a systematic mathematical calculation we have also to consider the energy dependence of ax, Eq. (35), in the 
expression ( 33 1 of the Fano asymmetry parameter associated with the first resonance, 



1 



1 

Tqx 



P'2 



2X 



E-Ex 



- F' 



(40) 



The function Xjqx is responsible for the asymmetry of the resonant contribution G^es, Eq. (31 1, that has a singularity 
at E = Eqx ■ The presence of a second resonance A' around A yields in 1 /qx a term singular at E = Eqx' ■ 

If the first resonance is very narrow and the second one broaden, <^ F^', the Fano asymmetry parameter l/qx 
varies slowly with the energy compared to the term in Gres singular at E = Eqx ■ In this case, the energy dependence 
of l/qx can be neglected around the resonance A and an energy-independent Fano asymmetry parameter can be 
defined. These results are in agreement with Ref. [19]. The resonant contribution to the conductance is then given as 
a Fano functiorP^ /(e) = |e + gFp/(e^ + 1) with a complex asymmetry parameter qp. For |l/q^F| ^ 1 this function 
has a quasi Breit-Wigner profile, while for II/^fI 3> 1 it becomes a symmetric dip, usually called antiresonance. 
The intermediate values |1/9f| — 1 correspond to a Fano function characterized by a maximum and a minimum 
approximatively equidistant to the axis e — 1 and we call this profile a "S-type" Fano line. For open quantum dots, 
the different Fano profiles can be associated with different types of interacting resonances. From Eq. (40 1 it follows 

that qx ^ Px- P'x foi' a second resonance much broader than the first one. The two vectors /3a and P'x, characterize the 
two considered resonances, and their scalar product is in principle nonzero. If the two resonant modes have different 
parities in the lateral direction, the vectors j3x and are approximatively orthogonal to each other and, in turn, the 
Fano asymmetry parameter has values from small to intermediate. In this case, we can speak about a weak interaction 
between the overlapping resonances. In contrast, for the same parity in the lateral direction, the vectors j3x and 
are approximatively parallel to each other and the Fano asymmetry parameter corresponds to a dip. In this case, the 
two overlapping resonances interact strongly. In Ref. [H] the antiresonances in the conductance through two identical 
quantum dots embedded in a wave guide were also related to strong interacting resonances with the same parity. 

In the case of two overlapping resonances with comparable widths, both the term in Gres singular at E = Eqx and 
the Fano asymmetry parameter ( 40 1 vary slowly with the energy and we can not predict the line shape around the 



resonances. This situation corresponds to a wide peak in the conductance. 

Summarizing all the above results and using the approximative expression of the total transmission around a 
resonance at the Fermi energy, Eq. (27), we obtain for the conductance 



GiVox + SV) ~ GresiEp - 5V; Vox) + G'^^siEr - SV; Vox) + GUEp - 6V; Vox) 



(41) 



where v^*-"""'""' = Vox is the potential energy in the dot region for which the resonance with the longest life-time 
{Tx < Tx') matches the Fermi energy. Based on the above relation, we identify the contribution of each resonance to the 
conductance and distinguish between weak and strong coupling regime of two overlapping resonances. The information 
about the strength of the coupling between the resonances A = {nx,ny) and A' = (n^,n^) is contained into the 
energy-dependent Fano asymmetry parameter, Eq. (33), and it determines the line shape of the resonant contribution 
Gres to the conductance. The other two components of the conductance, G'^^^ and G'^^, provide information about 
the interaction of the second resonance {n'^,n'y) with all other resonances of the system excepting the two already 
considered. A strong variation with the energy of the function G[,g in the energy domain of the resonance (n^,n^) 
indicates the presence of a third resonance (n", riy) around the Fermi energy. From the line shape of the resonant 
component G^g^ we can, in principle, get the information about how strong this third resonance interacts with the 
resonance (n'^^n'y). 



1. Weak interacting resonances 

In the weak interaction regime the two overlapping resonances are close in energy but they do not perturb each 
other significantly. Each of them contributes to the conductance as a quasi-isolated resonance and the line shape of 
the peak is given as a superposition of two Fano lines with a slight up to an intermediate asymmetry. This is the case 
of the second peak in Fig. [sja) and the first peak in Fig. |6]^a) , for which the different contributions to the conductance 
are analyzed in detail in Fig. flO] The two peaks correspond to the pair of resonances (1, 2) and (2, 1) and (2, 3) and 
(3,2). 

In both situations the overlapping resonances have at the origin a degenerate eigenstate of the isolated dot presented 
in Fig. [2]jb), with different symmetries in the x- and y-direction. The two quantum point contacts (V^i and Vh2) of 
the open dot create a strong coupling regime to the conducting leads and break the square symmetry of the isolated 
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FIG. 10: Conductance peaks associated with isolated resonances: (a)A = (1,2), A' — (2,1); (b)A = (2,3), A' = (3,2). Upper 
parts: Resonant parts of the conductance Gres and G^es and background part Gj,g; Lower part: Poles and position of the Fermi 
level in the complex energy plane. The potential energy in the dot region is constant, Vd ~ Vqa- 



dot. In turn, the degeneracy is lifted "when the quantum dot becomes open and, "with increasing the coupling strength, 
the degenerate energy level evolves into t'wo reson ance s that repulse each other in the complex energy pl ane. T his 
phenomenon is illustrated in the lowerpart of Fig. 10 The first case corresponds to the resonance trappin^^^ and 



the second one to the level repulsiotP^^^^. Due to the trapping, the resonance (1, 2) has a longer life-time and a resonant 
state almost localized inside the dot region, "while the resonance (2, 1) has a shorter life-time and sho'ws a significant 
probability distribution density in the region of the t'wo quantum point contacts. The state "with the nearest maximum 
to the aperture couples stronger to the contacts and yields a broader contribution to the conductance compared to the 
first resonance. Both of them are described by Fano lines "with a slight asymmetry corresponding to a "weak interaction 
bet"ween the overlapping resonances and "with the background. Figure 10 'b) shovifs the resonances (2,3) and (3,2), 
"which have comparable "widths and are "well separated in energy. From the probability distribution density of these 
t"wo modes. Fig. |6][a), it is evident that both of them can easy couple to the modes localized in the point contact 
regions, modes denoted by "a" and presented in Fig. |9] In this case, the resonance trapping is not favorable. The 
both resonant modes couple to the contacts via two modes of type " a" . The line shape of Gres corresponds to a "weak 
interaction bet"ween the resonances (2, 3) and (3, 2), but GJ.g^ indicates a stronger interaction of the second resonance 
"with the background. The total contribution to the conductance yields in each situation a peak "with a maximum 
about 2. 

In conclusion, the "weak coupling regime bet"ween overlapping resonances is characterized by probability distribution 
densities -within the dot region similar to the eigenstates of the isolated dot. According to this rule, the resonances 
(1,4) and (4,1) and (3,4) and (4,3) are also "weak coupled "with each other, but, as we "will see in the next section, 
in each case there is a strong coupling "with another neighbor resonance "which modifies the probability distribution 
density of the states (4,1) and (4,3), respectively. The last peak in Fig. IWa) corresponds to the resonances (2, 5) and 
(5, 2) that interact also "weakly and have a similar behavior to the pair (2, 3) and (3, 2). 



2. Strong interacting resonances. Hybrid modes 



The really ne"w physics of the scattering process can be seen in the case of a strong interaction bet"ween the 
overlapping resonances, phenomenon that does not occur in the case of an effective ID quantum dolj^^. This coupling 
regime is responsible for the thin or strong asymmetric peaks and dips in the conductance and for the resonant 
states "whose probability distribution densities differ strongly from the corresponding eigenstates of the isolated dot. 
Particularly for the SET geometry, Fig. [TJ the strong coupling of the quantum dot to the environment is al"ways 
accompanied by a strong scattering bet"ween the energy channels. This supplementary scattering determines the 
reordering process of the resonances in the complex energy plane, i. e., the interaction bet"ween overlapping resonances. 
The channel mixing infiuences especially the eigenstates "with the same symmetry in the lateral direction. Due to the 
favorable parity, these modes couple "with each other and generate ne"w resonant modes that can not be supported by 
the isolated dot. As seen in Figs. [5][6j and[7j there are t"wo categories of strong coupled resonances: The first ones 
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are the resonances that correspond to eigenstates with the same symmetry in the x— and y— direction and whose 
resonant states are hybrid modes, similar to the hybrid orbitals of the natural atoms. The resonances (1,3) and (3,1), 
(2,4) and (4,2), and (1,5) and (5,1) belong to this category. The second category includes resonances corresponding 
to eigenstates with the same symmetry only in the lateral direction (y-direction) like the pairs (4,1) and (3,3) and 
(4,3) and (5,1). We associate these resonances with a strong interaction because the probability distribution densities 
for the states (4,1) and (5,1) are drastically modified in comparison with the isolated case. 

The modes associated with strong interacting resonances yield dips or "S-type" Fano lines in the conductance, 
superposed on the top of broad peaks, as shown in Figs. [5|^ a), [6| ^a), and[7j^a). The overlapping resonance pairs (1, 3) 
and (3, 1) and (2, 4) and (4, 2) are analyzed in detail in Fig. |11[ The strong interaction of these resonances, reflected 
by the dips in Gres, has as an effect their strong repulsion in the complex energy plane. In turn, the resonant modes 
(1, 3) and (2, 4) become long-lived and are practically localized within the dot region, while the modes (3, 1) and (4, 2) 
couple stronger to the contacts and have shorter life-time. The electron distribution in the dot region favors in the 
first case, Fig. [ll|a), a typical resonance trapping, while in the second case. Fig. [lljb), this phenomenon is not so 
pronounced, but it is accompanied by a level repulsion on the real axis. The contribution of the second resonance to 
the conductance G'^.^^ is described by a broad peak and the background GJ, is almost constant. 
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FIG. 11: Conductance peaks associated with isolated resonances: (a) A = (1,3), A' = (3,1); (b) A = (2,4), A' = (4,2). Upper 
parts: Resonant parts of the conductance Gres and G'^^s ^md background part G^g; Lower part: Poles and position of the Fermi 
level in the complex energy plane. The potential energy in the dot region is constant, Vd ~ Vox 



Each pair of strong interacting resonances analyzed above corresponds to a degenerate energy level of the isolated 
dot and their probability distribution densities are practically linear combinations of the two eigenfunctions of the 
degenerate level. This property can be easy seen in Fig. [6f^a) for the resonances (2, 4) and (4, 2). We can speak in this 
case about hybrid resonant modes. The open quantum dot behaves like the oxygen atom in the water molecule: due 
to the interaction with the hydrogen atoms, the s and p orbitals of the oxygen are mixed to new hybrid orbitals so 
that the total energy of the molecule is minimal. Similar, the coupling of the quantum dot to the contacts by means of 
two quantum point contacts yields a supplementary scattering potential which allows for new resonant modes. They 
are not states which survive the coupling process to the contact^^, but rather new hybrid states, whose existence 
is directly connected to the presence of the strong coupling regime. These modes offer the possibility of engineering 
quantum systems with complex properties. Even in the case of a non-perfect square quantum dot the above results 
remain valid. A small difference between dx and dy yields instead of a degenerate level two very close eigenvalues. 
Essential for the strong interaction of the two corresponding resonances is the same parity of the resonant states on 
both directions and not the initial degeneracy. 

The resonances (1, 5) and (5, 1) interact also strongly. They determine in the conductance a very thin "S-type" Fano 
line superposed onto an extreme broad peak as shown in Fig. |12[ b). Their stronger repulsion in the complex energy 
plane compared to the precedent cases (Fig. 11 ) is determined by a supplemental strong interaction of the resonances 
(5, 1) and (4, 3), which have the same parity in the lateral direction. The maps of the probability distribution densities 
for the two modes in Fig. [Tf^a) confirm also the phenomenon of hybridization. The resonant modes (4, 3) and (5, 1) do 
not show such a high symmetry as the modes (2, 4) and (4, 2), but it is evident that they can be obtained as a linear 
combination of the eigenbmctions (4, 3) and (5, 1) of the isolated dot and the mode (4, 3) dominates this combination. 
The multiple interactions between neighbor resonances with the same symmetry in the lateral direction amplify the 
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FIG. 12: Conductance peaks associated with isolated resonances: (a) A — (1,4), A' — (3,3), A" = (4,1); (b) A = (1,5), 
A' = (3, 4), A" = (4, 3). Upper parts: Resonant parts of the conductance Gres, G'res, and G"es and background part G'^g\ Lower 
part: Poles and position of the Fermi level in the complex energy plane. The potential energy in the dot region is constant, 
Vd = Vox. 



phenomenon of resonance trapping. One resonance - in this case resonance (1,5) - decouples from the contacts and 
becomes extreme long-hved, while the other two become broaden and show a significant separation in energy. If we 
consider the quantum dot as an artificial atom it is easy to accept the hybridization as a natural process determined 
by the interaction with another system, but the existence of very narrow resonances supported by an open quantum 
dot seems to be a paradox: it is necessary to open a quantum dot, i. e., to allow for regions where the direct electron 
transfer between dot and contacts is possible, in order to obtain strongly localized states. Hence, long-lived modes 
of a quantum system can be obtained either in a quasi-isolated quantum dot or in a dot confined by shallow barriers 
engineered in such a way that the scattering channels are strongly mixed. The two types of localized modes have 
different fingerprints in the conductance: in the first case they yield quasi-symmetric maxima, while in the second 
case strong asymmetric Fano lines appear on top of broad peaks. 

The last sequence to be discussed corresponds to the resonances (1,4), (4, 1), and (3,3) in Fig. [6]ja). As can be 
seen from the probability distribution density maps in Fig. [6]^a), there are three interacting resonances with different 
coupling strengths; (1,4) and (4,1) interact weakly and yield two slight asymmetric maxima in the conductance, one 



of them quite thin and the other one broad, Fig. 12 ^a). In contrast, the resonances (4,1) and (3,3) interact strongly 



and the resonant contribution of (3,3) is a broad dip. The three interacting resonances (4,1), (3,3) and (1,4) are 
very interesting in view of the experiments presented in Ref. [5]. Their contribution to the conductance together 
with the next peak determined by the resonances (2,4) and (4,2), Fig. |6|^a), approximate qualitatively very well the 
conductance curve given in Ref. |6], Fig. 2(a), for the quantum dot in the Fano regime. Based on our resonance 
analysis we can conclude that the first thin peak in the measured conductance curve is superimposed on the top of a 
second broad peak and they correspond to two weak interacting resonances with different symmetries in the lateral 
direction. The next dip in the conductance reflects the presence of a resonance of type (n, n) that interacts strongly 
with only one of the neighbor resonances. The following "S-type" Fano line is again superimposed on a broad peak 
and indicates the presence of two strong interacting resonances with the same symmetry in the lateral and transport 
directions. For a quantitative analysis of the conductance we have to determine from the charge analysis within the 
dot region the value interval of Vd that corresponds to the number of electrons found experimentally^. 



V. CONCLUSIONS 



We have provided in this paper a systematic treatment of the conductance through a quantum dot strongly coupled 
to wide conducting leads via short quantum point contacts. The electronic transport through this type of dots is 
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essentially a scattering process in a low confining potential, which requires a direct solution of the two-dimensional 
Schrodinger equation with a nonseparable scattering potential. For this purpose, we have used a generalized scattering 
theory that allows for a complete description inside and outside the scattering area and is based on the R-matrix 
formalism. The resonances are determined as poles of the multidimensional scattering matrix, which contains the 
information about channel mixing due to the nonseparable scattering potential. The strong coupling of the quantum 
dot to the environment yields overlapping resonances, which show a significant interaction with each other in the case 
of a favorable parity of the corresponding resonant states. 

The conductance is determined as a function of the potential energy within the dot region, and every peak in 
conductance is associated with a resonance or a group of overlapping resonances. Based on the representation of 
the scattering matrix in terms of the R-matrix we provide for each peak an exact decomposition of the conductance 
in resonant terms associated with each of the overlapping resonances and a background. The decomposition is 
hierarchical, i. e., from the strongest to the broadest resonance, and allows for a deep understanding of the phenomena, 
which determine the transmission in the case of interacting resonances. The resonant states characterizing the open 
quantum dot in the Fano regime are presented in comparison with the eigenstates of the isolated dot. Every resonant 
state has a correspondent between these eigenstates and, we distinguish between slight and strong modified states 
due to the coupling with the environment. The last ones are called hybrid resonant modes, and they occur only in the 
case of a strong coupling regime of the quantum dot to the contacts, as an effect of the interaction between resonances 
with the same parity. The phenomenon of hybridization evidenced here for the quantum dots in the Fano regime 
of transport attests the molecule-like behavior of this system and opens the possibility to realize artificial molecules 
based on semiconductor nanostructures. 

The conductance through the quantum dot in the Fano regime of transport is also compared qualitatively to the 
experimental data reported in Ref. 
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Appendix A: Poles of the 5-matrix 

The starting point for our pole analysis is the expression of the non-constant part in the 5-matrix in terms of 
Ax^dx0l/{E-Ex) and fix: 

l + in = {l + iAx){l + inx). (Al) 



Using the definition of /3x, Eq. (22), we can immediately demonstrate that each determinant of the second order of 
Ax is zero, {Ax)ij{Ax)ip — {Ax)i^lAx)ij = 0, where each index l,p is a composite index (sn) with s = 1, 2, n > 1. 
Therefore, the matrix has the rank 1. On this basis we find that 

det [1 + iAx] = 1 + iTi-[Ax] . (A2) 

In order to demonstrate the above relation we consider a, M x M matrix A, M > 2, with Rank[A] = 1, calculate the 
determinant of 1 + A and take after that the limit M — oo. According to the definition, the determinant is given as 
a sum over all permutations n of the numbers {1, M} 

det[l + A] = ^SgnTTMiSlmi + Ai„i,)....{SMmM + ^MniAi)' (A3) 

with mi, m-M & {1, M}, mi ^ nij for i ^ j, i, j G {1, M} and sgn7rj\/ denotes the signature of the permutation 
TT and Sij is the Kronecker delta. After a direct calculation we obtain 



det [1 + A] = ^ sgn7rA/(5i„ii....(5 

+ X] SgnTTMfJlmi ••••<5A/-lmA/-i^ 



MrriA 
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Imi ^2m,2 ^Mm m 

fjii ^M—2mM -2 Am 



TIM 

+ All 



1+A 



^M-lAf-l ^Af-IM 
^AfAf-l ^AfAf 



-det[A]. 



Thus, det [1 + A] is given as 1 plus a sum of determinants of different order of A. But Rank[A] = 1 and consequently 
all the determinants of A up to the second order are zero. So that we find 



det [1 + A] = 1 + Tr[A] 



(A4) 



This result does not depend explicitly on the matrix dimension M so that we can generalize it for the case M — > cx) 
and obtain Eq. ( A2 ) . 



In the next step we calculate the adjugate matrix of 1 + A, i. e., 1 + A, in order to invert it: 

(i+Ari= ^ 



det[l + A]' 



(A5) 



For a given pair of indices ij the corresponding matrix element of 1 + A is calculated as the product between (—1)*+^ 
and the minor ij of (1 + A)'^ (the determinant of the matrix obtained from (1 + A)-^ by removing the row i and the 
column j, where T denotes the matrix transpose). Thus, for i = j we obtain 



(1+A),, = det[l + B,,], 



(A6) 



where B^^ is obtained from A"^ by removing the row i and the column i. The matrix 'Ba is a {M — 1) x {M — 1) 
matrix with the rank 1 and therefore 



A/-1 



(1 + A),, = 1 + Tr[B,,] ^l+Y. (^")« = 1 + Tr[A] 



In the case j — i + 1 wc find 



(1 + A)^^_^i = -det[Ii + B,,+i], 



(A7) 



(A8) 



where (M — 1) x (M — 1) matrix Ij is obtained from the unity matrix by changing 1 on the position ii with 0, 
(i < M) and 'Bu+i is the matrix obtained from A by removing the z-th row and the {i + l)-th column. This is a 
(M — 1) X {M — 1) matrix and has the rank 1. The matrix element ii of B^j+i is Aa-^i. Further we write explicitly 
the determinant as a sum over all permutation of the numbers {l,...,il/ — 1}, 

det[I,+B,,+i] = J2 

mi ^~ ^Imi inii 

TTJ .J - 1 

^{Si + lnii+i + fej+lmi+i )----(<5A/-lmA/-i +^A/-lmM-l)i 

where bji, with j,l = 1, Af — 1, means the matrix element jl of Bii+i. Replacing birm by Sirm + birm — ^inii allows us 
to express (1 + A)^^^-^ as det[l + Bii+i] minus the minor ii of I,; + Bii+i. The last two determinants can be calculated 
using Eq. ( A2 ) because the corresponding matrices are a sum of the unity matrix and a part of ^-matrix which has 
the property Rank[A] = 1. Thus, we find 



A/-1 

det[Ii + Bu+i] = 1+51 ^w' ^ 



j-l A/-1 
j = l j=i+l 



(A9) 



and after that, using Eq. (A4), 



(1 + A) 



ii+l 



-A,, 



(AlO) 
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Further we analyze the case j > i + 1. If we ehminate the i-th row and the j-th column in (1 + A)"^ we obtain a 
matrix which does not have any more elements of the type 1 + an on the main diagonal between the column i and 
j — 1. These elements are on the positions I — 11, I = i + 1, j — 2. Taking into account that we only need to calculate 
the determinant of this matrix we exchange the columns: ioi + l, j — 2-;->j/ — 1 and each of these j — i — I 
operations changes the determinant with -1. So that we can write 

+ = (-l)*+J"+J"-^-Met[I,_i + B,,]. (All) 

As described above the matrix B^j is obtained from A"^ by removing the i-th row and the j-th column and after that 
by exchanging the columns i ^ i + j — 2 ^ j — 1. This (M — 1) x (M— 1) matrix has also rank 1 and the matrix 



element j — Ij — 1 of B^j is A. 



1 ^ 



(1 



Similar to the previous case we can demonstrate here that 
A),, = -A^j. (A12) 



With the same procedure we can demonstrate that the above result remains also valid for j < i. 
If we put together the main results of this section, Eqs. (A2|, (A7), (AlO), and (A12), we find 



(l + A)-i = l- 



1 + Tr[A] ■ 

The above relation does not depend essentially on M, so that we can take the limit M 
for iA\. After that we obtain from Eq. (All that 



(A13) 



(i + in)-^ = (i + inx) 



oo and generalize Eq. (A13) 



(AM) 



Feeding this relation into the definition of 5-matrix, Eq. (15), we find Eq. (21 1. 
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